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ANALYSIS AND PREDICTION OF LONGITUDINAL STABILITY OF AIRPLANES 
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SUMMARY 

An analysis has been made of the longitudinal stability 
characteristics of 15 airplanes as determined in flight. 
In the correlation oj satisfactory and unsatisfactory char- 
acteristics with determined values^ the derivative ddjda 
that expresses the ratio of static-restoring moments to 
elevator-control moments was found to represent most 
nearly the stability characteristic appreciated by the pilots. 
In the derivative, de is the elevator angle and a is the angle 
of attack. ^ This derivative may be readily determined in 
flight or in the wind tunnel by measuring the elevator 
angles required for trim throughout the angle-of-attack 
range. It affords, moreover, a means of comparing 
airplane characteristics because the stick movement and 
the type of stick-force gradient are dependent on the 
magnitude of db^jda. 

The analysis was extended to study the effects of various 
design features on the observed stability characteristics. 
In this connection an expression was derived by means of 
which dbelda may be computed on the basis of general 
airplane dimensions for the propeUer-qff condition. 
Comparison of computed values, or comparable wind- 
tunnel determinations with engine-idling flight conditions, 
shows a powerful destahilizing efect of the idling propeller. 
An empirical expression based on the propeller dimensions 
brought the values obtained in the propeller-ofl tests and 
the computations into good agreement with the flight values. 
General power-on effects observed in the various airplanes 
tested are discussed. 

Design charts and data are included that show the effects 
on longitudinal stability of relative positions of wing and 
tail, fuselage size and location, engine nacelles, and hori- 
20 n tal-tail arra ngem e nts . 

Also included is a discussion of desirable numerical 
values of dde/da. For design a value of dhjda of 0.5 is 
suggested, 

INTRODUCTION 

A few years ago the National Advisory Committee for 
Aeronautics instituted a program for the study of the 
flying qualities of airplanes. The primary purpose of 
the research program was to determine quantitatively 
what constitutes satisfactory flying qualities and the 
stability and control requirements that an airplane can 
be expected to fulfill. Accordingly, the investigation 
has consisted mainly in determining these characteris- 
tics in flight tests of various airplanes. These airplanes 



were made available largely by the Army and more 
recently by private companies in accordance with 
requests by the Civil Aeronautics Authority. 

Drawings of the airplanes tested in flight are given in 
figure 1 and pertinent dimensions are listed in table I. 

Stability and control characteristics have now been 
determined for 15 airplanes of varied types. The data 
obtained have not only showTi what constitutes satis- 
factorv^ flying qualities but also show, by a proper 
analysis, how various design features influence the 
observed flying qualities. 

The present report presents the results of such an 
anatysis as regards the longitudinal-stability and 
control characteristics of the various airplanes tested. 

DEFINITION OF LONGITUDINAL STABILITY 

The longitudinal-stability characteristics of an air- 
plane are conveyed to the pilot as the variation of stick 
force and stick position with air speed or angle of attack. 
The pilot's opinion of the stability is therefore governed 
not only by the static pitcliing-momcnt characteristics 
dCJda of the airplane but also by the power of the ele- 
vators dCJdSe and their hinge-moment characteristics. 
Accordingly, a logical criterion for longitudinal stability 
is defined as the ratio of static-restoring moment to 
elevator-control moment dbjda because both the stick 
movement that the airplane experiences with a changing 
angle of attack and the stick-force gradient depend on 
this derivative. 

Specific values of this derivative desirable for the 
various types of airplane are still unestablislied. Flight 
tests indicate definitely, however, that the derivative 
should never be negative or of so low a value that a 
reversal of either the stick movcmetit or the stick force 
occurs. Preliminary analysis of the stick-force data 
shows that, for conventional tail arrangements, a value 
of dde/da of approximately 0.2 is required to insure 
stick-free stability because of normal elevator-floating 
tendencies with angle-of-attack changes. Stick move- 
ment is, of course, very nearly proportional to dSjda. 
On airplanes intended for high maneuverability, large 
values of dSjda have been found to be desirable pro- 
vided that they are obtainable without heavy control 
forces. For airplanes of this category, dSjda should 
be of such a value that at least a 4-incli movement of 
the top of the stick is required to change from a low 
to a high angle of attack in accelerated maneuvers. 
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Airplanes of appreciably lower stick travel have been 
found to be very susceptible of inadvertent stalling in 
maneuvers. 

No upper limit for db^jda has been indicated except 
by considerations of maximum control requirements, 
such as in landing with flaps extended. It may also be 
true that for airplanes of the transport and the heavy- 
bomber type a sufBcient hinge-moment reduction cannot 
be obtained to allow appreciably greater values of 
dbeldoc than arc representative at the present time. For 
reasons such as these, specific values cannot be recom- 
mended with assurance although all flight data indicate 
that large values of dde/da are in themselves desirable 
provided that they are not obtained at the expense of 
acceptable control-force characteristics. 

Additional discussion concerning numerical values of 
ddjda is included later in this report in connection with 
suggested design practice. 

STABILITY EQUATION 

Preliminary^ attempts to correlate the stability char- 
acteristics determined in flight for the engine-idhng 
condition with those computed on a basis of wing-and- 
tail theory alone indicated the necessity of considering 
the destabilizing effects of the fuselage, the nacelles, and 
the idhng propeller. In the initial development of the 
stability equation the effect of the idling propellers is 
ignored. 

Consideration of the forces and tlie moments acting 
on an airplane without propeller leads to the following 
expression for dbjda: 



d8e__ 
da 



where 



2o dai 



da 

Suf 



K, 

If 

N 

ij 



elevator effectiveness factor 

KdC^JdaJ 

rate of change of dowmwash over horizontal tail 

with angle of attack 
wing area, including section through fuselage and 

ailerons, square feet 
horizontal distance from aerodynamic center of 

wing to airplane center of gravity, feet 

slope of wing Hft-coefficient curve, per radian 

fuselage and engine nacelle-moment factor 
maximum fuselage width, feet 
over-all fuselage length, feet 
number of engine nacelles in a multiengine air- 
plane 

maximum engine nacelle width, feet 

over-all length of engine nacelle, estimated to be 

streamline body, feet 
ratio of dynamic pressure over horizontal tail to 

free-stream dynamic pressure (0.9) 
horizontal distance from airplane center of gravity 

to elevator hinge line, feet 



St total horizontal tail area including section through 

fuselage, square feet 
dCs^ slope of curve of normal-force coefficient for hori- 
dai zontal taU, per radian 
and 

A aspect ratio 

r factor in expression for slope of normal-force curve 
for tail surfaces with end plates 
The defined symbols that occur in equation (1), the 
figures that are to be used for their evaluation, and the 
references from which the data for the figures were 
taken are listed in table II. 

Table XL— REFERENCES FOR EVALUATION OF 
SYMBOLS USED IN EXPRESSION FOR dBJda 
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A7 
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' For horizontal tails with vertical tails mounted as end plates, see figure 9. 




Figure 2.— Dimensions used In stability calculations. Distances measured at anple 

of zero lift. 
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iMost of the quantities in equation (1) may be com- 
puted by the use of the figures listed in table II. Some 
explanation of the use of the charts in figures 4, 5, and 
6 is, however, necessary for the evaluation of dtlda. 

The procedure to be followed in the determination of 
dtlda may be summarized in the follo^ving steps: 




./ .a .J , 

Elevofor area back of hinge line Se 
Horizon f of foil area 



Figure 3.— Variation of elevator effectiveness factor r with ratio of elevator area 
back of hinge line to horizontal tail area. Data taken from reference 1. 



(1) Determine which group of the four groups given 
in figure 4 most nearly approximates the design taper 
ratio and aspect ratio. This step is not critical because 
corrections must later be applied for diflerences in aspect 
ratio and taper ratio. 

(2) In the selected group, using values of mo, x, and 
Xx as defined in figure 2, with linear interpolations for x 
and Xt wherever necessary, determine rfe/rfa. 

(3) In figure 5 (a) plot d^lda as found in step (2) 
against the group taper ratio and from this point follow 
the contour line to the design taper ratio and read of! a 
new value of dtjda. 

(1) In figure 5 (b) determine the correction to dtjda 
for the difference in the aspect ratio between the group 
value and the design value at the design taper ratio. 
Add this correction ^{dtlda) to the value of rfc/rfa de- 
termined in step (3) if the design value of the aspect 
ratio is low^er than the group value; subtract this cor- 
rection if the design value of the aspect ratio is higher 
than the group value. 

(5) From figure 6 determine the correction factor to 
be multiplied by the total value of dtjda to convert it 
from the absolute value at the center line to an average 
value over the horizontal tail span. 

With the following additional dimensions an example 
of the computation of dtjda for airplane 1 (fig. 1) is indi- 
cated in figures 4, 5, and 6. 
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Taper rafio _ Tail span 

(a) Corrections to downwash for variations in taper ratio. W(ng span 

(b) Corrections to downwash for variations in aspect ratio. 

FinuRE 5 -Corrections to downwash for variations in wing aspect and taper ratios. Figure C.-Correction factor to convert downwash at center of horizontal tail to 

Data taken from reference 2. average downwash over tail. Data taken from reference 2. 




Figure 7.— Variation of dCs^fdat and dCJda^ with aspect ratio for conventional sections. Data taken from reference l. 
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Horizontal distance from wing root trailing edge to 
hinge line, (6/2) feet 

Horizontal distance from wing root }4-chord 
hinge line, (x)(6/2) feet 

Vertical distance from wing root trailing edge 
line, (mo)(b/2) feet___ 

Tlien 

Wo=;;^ = 0. 084 



elevator 

30.5 

point to elevator 

51.9 

to elevator hinge 
6. 3 



74. 5 
30.5 



74. 5 



= 0.41 



3-= 



51. 9 
74.5 



=0. 696 



3.2 
2.8 
2.4 
20 
IS 
1.2 
£ 



0 20 40 60 

Wing section 1/4 -chord post Hon, 
percent body lengtti 

Figure 8.— Variation of fuselage and engine nacelle moment factor Ks with relative 

location of wing and body. Data taken from reference 3. 

These values determine point A in figure 4 (step 2). 
If point A is transposed to figure 5 (a) and the contour 
lines are followed, as indicated by the dashed line, point 
B is located (step 3). The increment BC, transferred 
from figure 5 (b) to figure 5 (a), establishes a value of 
dtlda of 0.56 (step 4). 

Interpolation is necessary in figure 6 to determine 
the correction factor of 0.95, which, when multiplied by 
the value of dtlda of 0.56, yields 0.53 (step 5). 

The discontinuities evident in the curves of figure 4 
arise from the fact that the empirical expression for 
dtlda. used in constructing the curves contained an abso- 
lute value of a for which a mean value was assumed. 
The discontinuities then occur for tail locations that 
bring the horizontal tail on the wing w^ake center at the 
assumed angle. 

Actually, because the wake has a finite thickness 
and because other locations near those corresponding 



to the discontinuities will similarly afl'cct the tail, the 
curves should be faired in some way. This fairing has 
not been attempted because the results would be mis- 
leading; if the tail traverses the wing wake at some por- 
tion of the angle-of-attack range, it will be operating 
at different values of dtlda before and after passing 
through the wake. Strictly, then, no single faired 
value of dtlda is applicable throughout the speed range 
for these tail locations. 




■'o ./ .a 

End-plafe t^eighf 
Horizontal tail span ' 6^ 
Figure 9.— Correction for dCstfdat of horizontal tail with end plsites. 

dtxt 




.4 .6 
d6g/dof (compufed) 

F[GUEB 10-— ComparlsoTi of experimental and computed values of stability, propeller 

off. 

The charts of figure 4 also indicate qualitatively the 
effects on dh^lda of wing dimensions and wing and tail 
arrangements. Thus the increased stability of a high- 
wing-monoplane arrangement over that of a low-w4ng 
monoplane may be partly attributed to the relative 
vertical location of the wing and the horizontal tail. 
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A comparison of propcllcr-off wind-tunnel results 
with computed values of db^jda is given in figure 10. 
The feasibility of predicting propcllcr-off stability 
characteristics by means of the expression given is indi- 
cated by the relatively small deviations of the data from 
the line of perfect agreement. 

Values of dbjda for wind-tunnel models for which 
onl}^ curves of dCm/da are available may be calculated 
from the following expression: 

dde_ Su>Cu, dCm 

dot 7 o Qt dCf^* da 

ThOt— — 

2o da, 

where is the mean aerodj^namic chord of the wing. 
CORRELATION WITH FLIGHT RESULTS 
A comparison of the computed values of dSjda with 
the values obtained in gliding flight with idling pro- 
pellers is showTX in figure 11. A destabilizing moment 
is clearly evident in the lack of agreement. 




O .^ .4 .6 ,8 i.O 

ddg/dcC (computed) 

FiGXTRE IK— Comparison of flight stability as determined with idling propellers with 

computed sttibility, propeller olT. 



It was found that the effect of the idling propeller 
could be represented by a term composed of the pro- 
peller dimensions with a constant empirical coefficient: 

with previously undefined symbols listed below: 
Kp empirical propeller coefficient (0.65) 
Np number of propellers 
D diameter of propeller 

horizontal distance from center of gravity to propeller 
plane, feet 
D and Ip are given in figure 2. 

The value of ddjda for gliding flight with idling 

propellers then becomes: 

r ^ JC. 



A comparison of the values of d8e/da as computed 
from this expression with those determined in flight is 
given in figure 12, 
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dS^jdoC (computed) 

Figure 12.— Comparison of computed stability with flight stabiliry as determined 
with idling propellers. 

LIMITATIONS OF METHOD 

Some deviation from perfect agreement may be noted 
hi figure 12. Several possible explanations for these 
irregularities arc: 

(1) The effect of the vertical center-of-gravity posi- 
tion has been neglected. This procedure makes the 
method somewhat conservative for high- wing mono- 
planes, particularly at high angles of attack. For low- 
wing and mid wing monoplanes the vertical center of 
gravity is usually sufficientl}' close to the aerodynamic 
center to be neglected, at least for attitude changes of 
the magnitude cxpericncod in power-off flight. For 
power-on flight, however, the vertical center-of-gravity 
position has a marked effect because of the large atti- 
tude changes experienced in going from low to high 
angles of attack. 

(2) As calculated, the value of dtjda is substantially 
independent of Cx. Under certain conditions, usually 
near the stall, a local flow brcakdowTi may make the 
computed results inapplicable. In many airplanes this 
breakdown results in values of dh^lda near the stall 
that are much greater than the values calculated or 
existing over the rest of the angle-of-attack range. 

(3) All of the effects, of the idling propeller cannot 
be adequately expressed by the simplified empirical 
term used. In this connection it should be appreciated 
that the value of the coefficient for the propeller- 
moment term was obtained for conventional tractor 
arrangements and may not be applicable for pusher or 
unconventional arrangements. 
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(4) Because few data are available for the detemii- 
nation of Kf at the more rearward positions of the wing 
with respect to the body, these values are subject to 
modifications. In addition, no allowance was made for 
the cross-sectional shape or plan form of the fuselage 
and the nacelles in the calculation of their moments. 

(5) The flight values ot ddjda, although substantially 
constant over the normal-flight range, w^ere subject to 
some individual variations. In the present analysis 
minimum values were always taken. 

POWER EFFECTS 

For all the airplanes tested, the power-on and the 
engine-idling values of d8e/da were substantially the 
same at low angles of attack. Without exception, low- 
ing monoplanes progressively lost stability with power 
on as the angle of attack increased. A large part of 
this effect was merely an increase in the power of the 
elevators as a result of the increase in qtfqo that occurred. 
This conclusion is indicated by the fact that the stability 
lost was a function of the original power-off stability. 
Low-^\ing monoplanes, however, that were neutrally 
stable with engine idling or with power on at high speeds 
also lost stability, a representative amount for high 
angles of attack being approximately dde/da=^— 0.2, 

On the basis of present knowledge, placing the center 
of gravity below the aerodynamic center appears to be 
the most direct method of limiting the stability loss due 
to power. This effect can be easily computed from 
considerations of the power-on attitude changes. High- 
wing monoplanes tested in flight showed relatively little 
stability loss due to power; midwing monoplanes were 
generally intermediate in this respect, 

SUGGESTED DESIGN PRACTICE 

Observations indicate that a value of dde/da of 
approximately 0.5 should be set as a lower limit of de- 
sirable stability in the gliding condition. Experience 
with several pursuit-type airplanes has shown that, 
with conventional ratios of stick to elevator movement, 
this value of propeller-idling stability results in approxi- 
mately 4 inches of stick movement in power-on ma- 
neuvers. With the center of gravitj^ well forward, this 
value of dhjda may be obtained with a small tail 
volume. This procedure is not recommended, however, 
because it results in undue sensitivity to changes in the 
ccnter-of-gravity location and makes more critical the 
control requirements for three-point landings. It is 
better to use a large tail volume with normal center-of- 
gravity locations. 

For large airplanes not required to maneuver, in 
which either visual or instrument references are always 
available to the pilot, this stability limit does not appear 
as essential. It is always desirable, however, that the 
stability be sufEciently great that a reversal of stick 
travel or stick force does not occur. 

The value of dhlda of approximately 0,5 is dictated 
not alone by the fact that large values are in themselves 
desirable; in addition, the reduction in stabihty to be 
anticipated with power on at low speeds should not 
diminish dhjda to undesirably low values. 



The methods of predicting longitudinal stability 
described in this report do not include the effects of 
wing flaps. For all the airplanes tested, however, the 
flap-down engine-idling conditions showed improved 
stability over that obtained with the flap up although 
with the flap down the destabilizing effect of power was 
in some cases greater than with the flap up. 

The requirements for stick-free stability are met 
when the airplane has only one trim speed for a given 
trim tab setting, above which speed push forces are 
required and below which, pull forces are required. A 
positive gradient of stick force against angle of attack 
of a magnitude sufficient to minimize the effects of con- 
trol friction sliould exist throughout the speed range. 
Preliminary investigation of the flight data indicates 
that a value of dddda of about 0.2 must be exceeded to 
obtain stick-free stability. 

Other considerations, chiefly the ability to lower the 
tail to the three-point attitude for landing with flaps 
down, determine the upper limit for dde/da, although it 
should be apparent that, where large values of dSjda 
are obtained by increased tail volume with a conven- 
tional amount of elevator control, no difficulty should 
occur. 

CONCLUDING REMARKS 

As a result of the foregoing analysis of flight data of 
longitudinal stability, the following conclusions may be 
stated: 

1. For many purposes the criterion of longitudinal 
stability is logically taken as dbjda, the rate of change 
of elevator angle with angle of attack, where 5e is the 
elevator angle and a is the angle of attack. 

2. From only a knowledge of the basic airplane di- 
mensions, dbjda may be predicted for both propeller- 
off and propeller-idling conditions for conventional 
airplane designs. 

3. The empirical factor that expresses the destabiliz- 
ing effect of the idling propeller may be used to correct 
propeller-off wind-tunnel data. 

4. For all airplanes tested the stability with the pro- 
pefler idling and with power on at cruising and high 
speed were very nearly tlie same. Low-\\ing mono- 
planes progressively lost stability as the angle of attack 
increased to the stall and high-wing monoplanes tended 
to retain their propeller-idling stability. 

5. For design a value of dbjda, with the propeller 
idling, of 0.5 is suggested to insure power-on stability 
and to give adequate stick movement for airplanes in- 
tended for high maneuverability. 

6. The suggested values of dbeldoc should be obtained 
by an adjustment of the effective tail volume with the 
normal center-of-gravity positions and the conventional 
amounts of elevator control. 



Langley Memorial Aeronautical Laboratory, 
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Table L— CHARACTERISTICS OF AIRPLANES TESTED 



Air. 
plane 
num- 
ber 


Wing 
area, 

(sq ft) 


Wing 
span, 

(ft) 


Distance 
from 
a c to 

"(ft)" 


Aspect 
ratio of 
wing, 
A 


Taper 
ratio of 
wing, 
X 


norizon- 
tal tall 
area, 

s, 

(sqft) 


Horizon, 
tal tail 
span, 
6i 
(ft) 


Aspect 
ratio of 

tail, 

A, 


St 


JJlSttiTlCe 

from 

C. g. to 

elevator 
hinge, 

(ft) 


1 


2780 


149 


-0. 45 


8.0 


4. 33 


505 


45.0 


4.0 


0.36 


49.5 


2 


1420 


104 


-.45 


7.6 


2.66 


254 


33.8 


4.5 


.37 


38.8 


'S 


965 


89.5 


.77 


8.4 


2.78 


203 


26.7 


3.5 


.30 


37.9 


3 


965 


89.5 


-.29 


8.4 


2.78 


203 


26.7 


3.5 


.30 


38.9 


3 


965 


89.5 


-.66 


8,4 


2. 78 


203 


20.7 


3.5 


.30 


39.4 


4 ,„ 


678 


70.5 


-.21 


7.3 


2. 78 


98 


20.2 


4.2 


.33 


28.5 


6 


545 


65.5 


-.58 


7.9 


3, 45 


116 


21.7 


4.1 


.34 


28.5 


6 


352 


49.5 


-.37 


7.0 


3.00 


87 


18.3 


3.9 


.39 


23.4 


7 


224 


36 


-.26 


5.8 


(-) 


41 


11.5 


3.2 


.33 


17.2 


8 


220 


36 


-.14 


5.9 


(-) 


41 


11.5 


3.2 


.33 


15.9 


g 


305 


42 


0 


5.8 


1.50 


47 


13.3 


3.8 


.40 


21.7 


10 


258 


39 


-.28 


' 5.9 


1.50 


60 


14.8 


3.6 


.46 


17.1 


u 


236 


37 


-.31 


5.9 


2.50 


48 


12.8 


3.4 


.3.^ 


17.9 


12 


248 


42 


.07 


7.1 


2.00 


48 


13.0 


3.5 


.39 


16.6 


13 


155 


34 


-.08 


7.5 


1.00 


28 


9.3 


3.1 


.39 


13.7 


13 


155 


34 


-.22 


7.5 


1.00 


28 


9.3 


3.1 


.39 


13.6 


13 


155 


34 


-.40 


7.5 


1.00 


28 


9.3 


3.1 


.39 


13.4 




1G9 


36 


-.15 


7.7 


1.00 


25 


10.2 


4.1 


.46 


13.6 


15.-_ 


180 


36 


-. 13 


7.2 


1.00 


26 


10.0 


3.9 


.42 


14. 1 



Fuselage 
length, 

Lf 

(ft) 



87.6 
68.0 
56.0 
56.0 
56.0 
44.7 
43.7 
36.3 
25.6 
25.0 
34.0 
28.2 
29.0 
27.6 
22.0 
22.0 
22.0 
18.5 
19.8 



Maxi- 
mum 

fuselage 
width, 

Wf 

(ft) 



10.4 
8.5 
0.3 
0.3 
0.3 
3.7 
5.5 
5.4 
4.5 
4.5 
3.8 
4.0 
4.4 
3.9 
3.5 
3.5 
3.5 
3.3 
3.5 



Equiva- 
lent na- 
celle 
le^th, 

(ft) 



25.0 
20.0 
15.4 
15.4 
15.4 
11.2 
16.0 
10.7 



Maxi- 
mum 
nacelle 
width, 

Wn 

(ft) 



6.4 
4.3 
4.3 
4.3 
4.3 
4.0 
4.0 
3.9 



Propeller 
diameter, 
D 
(ft) 



11.5 
11.5 
11.5 
11.5 
11.5 
11.5 
10.5 
8.8 
10.5 
10.5 
11.0 
9.5 
10.0 
9.0 
6.0 
6.0 
6.0 
5.8 
6.5 



Distance 
from 

c. g. to 
propeller 

plane, 

(ft) 



13.4 
13.9 
11.6 
12.0 
13.1 
7.4 
8.0 
C.4 
6.5 
7.1 
8.9 
7.8 
7.0 
6.4 
6.2 
6.3 
6.5 
4.7 
5.3 



■ EUiptlcal. 



